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^ ■ ABSTRACT 

CZ5 . 

Aims. We present observations and analysis of time-series spectroscopy and photometry of the pulsating subdwarf 
B star PG 1219+534 (KY UMa). Subdwarf B stars are blue horizontal branch stars which have shed most of their 
hydrogen envelopes. Pulsating subdwarf B stars allow a probe into this interesting phase of evolution. 
Methods. Low resolution spectra were obtained at the Nordic Optical Telescope and Kitt Peak National Observatory, 
and photometric observations were obtained at MDM and Baker observatories in 2006. We extracted radial velocity 
and equivalent width variations from several Balmer and He I lines in individual spectra. The pulsation frequencies 
were separated via phase binning to detect line-profile variations in Balmer and helium lines, which were subsequently 
matched to atmospheric models to infer effective temperature and gravity changes throughout the pulsation cycle. 
Results. From the photometry we recovered the four previously observed frequencies and detected a new fifth frequency. 
From the spectra we directly measured radial velocity and equivalent width variations for the four main frequencies 
and from atmospheric models we successfully inferred temperature and gravity changes for these four frequencies. We 
compared amplitude ratios and phase differences of these quantities and searched for outliers which could be identified 
as high-degree modes. These are the first such measurements for a "normal" amplitude pulsating subdwarf B star, 
indicating that spectroscopic studies can benefit the majority of pulsating subdwarf B stars. 



Key words, stars: subdwarfs - stars: variable: general, subdwarf stars - stars: individual: PG 1219+534 - Techniques: 
spectroscopic, radial velocities 

1. Introduction of the same mass (to within a few perce nt) every time 

(iD'Cruz et al.l . ll996UHan et al.l . l2002l . l2003h . Pulsating sdB 
Subdwarf B (sdB) stars are the field counterparts to ex- stars come in tw0 varie ties. The short perio d (90 to 600 sec- 
tended horizontal branch stars observed in globular clus- onds ) pu l sa tors were discovered in 1997 (|Kilkennv et all 
ters. Their masses are about 0.5 M Q with thin (<10 M Q ), [199^ and were named EC 14026 sta rs after that pro- 
inert hydrogen envelopes and temperatures from 22 000 to toty pe. Their official designation is V361 Hya stars, the 
40 000 K (Heber et al. 1984; Saffer et al. 1994), making General Catalog o f Var iables Stars lists them as RPHS 
them very blue. Their origin remains a mystery as mass (|Kazarovets et all . [20031) . but they are typically referred 
loss on the first-ascent giant branch must remove all but to as sdBV stars xheir usual amplitudes are near 1% of 
< 10 M Q of the hydrogen envelope, yet produce a core their mean brightness and s tars have be en detected with as 



few as one (HS 1824 +5745 dReed et alll2006D ) or as many 



Send offprint requests to: M. D. Reed as 55 (P G 1605+072 (Kilkenny et all . 11999)) pulsation fro 

* Based on observations from the Nordic Optical Telescope, quencies (jReed et ail . l2007al) . Long period (4 5 minutes to 

Kitt Peak National Observatory, and MDM Observatory. 2 hours) pulsators were discovered in 2003 (jGreen et all 

** Visiting Astronomer, Kitt Peak National Observatory, [2003; Ree d et al.l . l2004a[) and are officially designated as 

National Optical Astronomy Observatory, which is operated V1093 Her stars but commonly named PG 1716 stars after 

by the Association of Universities for Research in Astronomy that prototype, PG 1716+426. Their amplitudes are typi- 

(AURA) under cooperative agreement with the National Science cally <0 .l% and they tend to be multimode pulsators. For 
Foundation. 



2 



M. D. Reed et al.: Time-series Spectroscopy of PG 1219+534 



this work, our interest is the sdBV class of pulsators as 
their periods are short, so many pulsation cycles can be 
observed during each run from a single site. The short pe- 
riods also produce spectroscopic variations (radial velocity, 
equivalent width, effective temperature and gravity) at an 
observable level for intermediate resolution spectrographs 
and 4 m class telescopes. 

Pulsating sdB stars potentially allow the opportu- 
nity to discern their interior structure using asteroseis- 
mology to obtain estimates of total mass, luminosity, en- 
velope mass, radiative levitation, gravitational diffusion, 
and helium fusion cross sections. To do so, the pulsa- 
tion modes must first be identifie d with their associ- 
ated spher ical harmonics. In 1999, 10 'Toole et ail (|2000t 
2002, 2003) began applying time-series spectroscopic tech- 
niques to attempt mode identification in sdBV stars. 
Whereas photometry measures brightness changes largely 
caused by temperature variations, spectroscopy can re- 
veal information regarding the pulsation velocities and 
separate the temperature and gravity components. An 
excellent example of th ese techniques was shown in 
iTelting fc 0stensenl (|2004f ) where they determined that the 
pulsational degree of the highest-amplitude frequency of 
PG 1325+101 (QQ Vir) was consistent with I = 0. Other 
sdBV stars studied using time-series spectroscopy include 
PB 8783 (EO Cet) and K PD 2109+4401 (V2203 Cyg) 
Jefferev fc PollaccoL 120001), PG 1605+072 (V338 Ser ) 
O'Toole et all . l2000Ll2002f 2003. 2005; Wool f et al.l.l2002f) . 



and Balloon 090100001 (|Telting fc 0stensenl . [2006) us- 
ing low-resolution spectroscopy on intermediate-sized tele- 
scopes. PG 1219+534, PG 1605+0 72, and PG 16 1 3+42 6 
were also observed using FUSE (|Kuassivi et all I2005T) . 
Only for the high-amplitude pulsators (PG 1605+072, 
PG 1325+101, and Balloon 090100001) were pulsations de- 
tected. However, such pulsators provide their own prob- 
lems in that the high-amplitude pulsations affect the other 
frequencies and therefore req uire considerable prewhiten - 
ing (variation rem oval); e.g. ITelting fc Ostensen] (|2006h : 
iTillich et alj (|2007f ). Unfortunately, because of insufficient 
signal, temporal resolution, or run length, none of these 
studies have provided conclusive mode identifications which 
would be useful for constraining models. 

This paper reports our time-series spectroscopic and 
photometric observations of PG 1219+534 (also KY UMa 
but hereaft er PG 1219) PG 1219 was discovered to be a 
pulsator by iKoen et ail (|1999f ) who detected four indepen- 
dent frequencies. The spacings of the frequencies are such 
that the entire pulsation spectrum can be resolved in just 
a few hours, yet the frequency density is too large for the 
pulsations to be of the same modal degree I. A photometry- 
based mode identification has been published which at- 
tribut es the four frequencies to three differing modal de- 
grees ()Charpinet et all l2005f ) which can be used as a guide 
for our mode identifications. While the pulsation ampli- 
tudes are low (2-8 mma in photometry), they should not 
interfere with each other, eliminating the need to prewhiten 
some frequencies to reveal others. This makes PG 1219 an 
excellent target for mode identification studies. Our pho- 
tometric observations of PG 1219 began in 2003 as part 
of our effort to resolve th e pulsation spectra (Fourier tran s- 
form; FT) o f sdBV s tars dReed et~aill2004bl [21M l2007allbl : 
IZhou et all I2006D ; lHarms et all (|2006f ) reported the re- 
sults of 2003 - 2005 photometric observations of PG 1219. 
During that time, the pulsation frequencies remained con- 



sistent with those of IKoen et alj (|1999[ ) and occasionally a 
new, low-amplitude frequency would appear for one or two 
nights. 

In §2 we describe our simultaneous spectroscopic and 
photometric observations during 2006 and examine the fre- 
quency content in §3. In §4 we discuss the results of our ob- 
servations and in §5 we provide conclusions from our work. 
Further detailed modeling to understand our observational 
results in terms of asteroseismology will be presented in a 
separate paper. 

2. Observations 

2.1. Photometry 

We obtained photometry at MDM and Baker observato- 
ries to support our spectroscopic observations, the details 
of which are provided in Table [1] Baker observatory (BO) 
is equipped with a 0.4 m telescope and a Roper Scientific 
RS1340b CCD photometer. At MDM we used the 1.3 m 
McGraw Hill telescope with an Apogee Instruments U47 
CCD. For all of our CCD measurements, we binned the 
CCDs 2x2 pixels providing a dead time of 1 second and 
all used a red cut-off filter (BG38 or BG40), so the effective 
bandpass covers the B and V filters and is essentially that of 
a blue-sensitive photomultiplier tube. Accurate timing was 
accomplished via Network Time Protocol (NTP) connec- 
tions to at least one stratum 1 and two stratum 2 servers. 
NTP timings were accurate to better than 0.02 s. 

Our photometric data were reduced using standard 
IRAF packages for image reduction, including bias sub- 
traction, dark current and flat field corrections. Intensities 
were extracted using IRAF aperture photometry with ex- 
tinction and cloud corrections using the normalized in- 
tensities of several field stars. As sdB stars are substan- 
tially hotter than typical field stars, differential light curves 
are not flat due to differential atmospheric and colour ex- 
tinctions. A low-order polynomial was fit to remove these 
trends from the data on a night-by-night basis. Finally, the 
lightcurves are normalized by their average flux and cen- 
tered around zero so the reported differential intensities are 
AI = (1/(1)) — 1. Amplitudes are given as milli- modulation 
amplitudes (mma) with an amplitude of 10 mma corre- 
sponding to 1.0% or 9.2 millimagnitudes. A portion of data 
from MDM is shown in Fig. [T] There is obvious beating that 
occurs on multiple time scales which indicates the multi- 
periodic nature of the pulsations. 

2.2. Spectroscopy 

During 7 nights in April 2006 we obtained nearly 5200 
time-resolved low-resolution spectra of PG 1219, with the 
specifics of each run provided in Tabled We obtained 2085 
spectra at the Kitt Peak National Observatory's (KPNO) 
4 m Mayall telescope using the RC spectrograph in long-slit 
mode. Our KPNO instrumental setup was: grating KPC- 
007 with the T2KB CCD which sampled approximately 
3000-5000 A, a FWHM spectral resolution of about 4 A, 
and a dispersion of 1.44 A/pixel. We subframed the CCD 
along the slit, reading out a section 2080 x 37 including 
a 41 column overscan region. Integration times were 10 s 
with rather long dead-times of 25 s for a total cycle time 
of ~ 35 s. Binning the CCD provided only negligible im- 
provement, so we did not. We obtained helium-neon- argon 
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Table 1. 2006 Photometric observations of PG 1219 simultaneous or near to the spectroscopic observations. 
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Run 


Date 


Start Time 


Length 


Number of 


Observat. 


Name 


(UT) 


(UT) 


(hours) 


images 




mdm040806 


8 Apr. 


3:10:30 


9.3 


3900 


MDM 


mdm040906 


9 Apr. 


3:21:30 


4.4 


1916 


MDM 


mdm041006 


10 Apr. 


2:24:52 


2.8 


521 


MDM 


mdm041106 


11 Apr. 


4:23:30 


3.5 


1132 


MDM 


bak041806 


18 Apr. 


5:09:30 


4.9 


1583 


BO 


bak041906 


19 Apr. 


2:23:30 


7.7 


2056 


BO 


bak042006 


20 Apr. 


2:20:00 


6.6 


1815 


BO 


bak042206 


22 Apr. 


2:25:00 


7.4 


2030 


BO 
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Fig. 1. Photometric lightcurves for consecutive 1 hour seg- 
ments of PG 1219. Solid line is a four-frequency fit to the 
points. 



Table 2. Spectroscopic Observations of PG 1219. 



Date 
2006 


Start Time 
(UT) 


Length 
(hours) 


Number of 
spectra 


Observat. 


8 Apr. 


3:23 


8.9 


810 


KPNO 


9 Apr. 


2:47 


9.5 


850 


KPNO 


11 Apr. 


7:34 


4.7 


425 


KPNO 


17 Apr. 


20:43 


8.3 


757 


NOT 


18 Apr. 


20:32 


8.9 


760 


NOT 


19 Apr. 


20:29 


8.7 


826 


NOT 


20 Apr. 


20:31 


8.0 


756 


NOT 



arc- line spectra every 80-150 spectra aperiodically and ob- 
served the flux-standard star Feige 66 three times per night. 
Halogen lamp flat fields were obtained nightly for removing 
pixel-to-pixel variations on the CCD. 

We obtained 3099 spectra at the Nordic Optical 
Telescope (NOT) using ALFOSC in long-slit mode. Our NOT 
setup used grism #16 and CCD #8, which sampled approx- 
imately 3500-5050 A, a resolution of about 3 A FWHM 
at a dispersion of 0.77 A/pixel. The grism and slit were 



set up such that they were aligned with the rows of the 
CCD to shorten the dead times. Using 25 s integrations, 
we achieved a cycle time of ~31 s, so while the NOT in- 
tegrated longer than KPNO, the cycle time was similar. A 
few re-acquisitions were done per night to correct the slit 
angle with the parallactic angle. Thorium-argon and helium 
arc-line spectra were obtained every 60-80 spectra. Halogen 
lamp images were obtained for flat-fielding. 

The spectra were bias and dark-current corrected us- 
ing overscan regions and flat field corrected using standard 
tasks within iraf. Two bad columns in the NOT data 
were corrected by linear interpolation of pixels in adja- 
cent columns. As the spectra contains low signal-to-noise 
(S/N) regions in the UV, and the KPNO CCD had a sharp 
sensitivity decline redward of 5030 A, the spectra were 
trimmed to 3710-5130 A for the KPNO data and 3480- 
5010 A for the NOT data. One-dimensional spectra were 
optimally extracted after subtracting a fit to the sky back- 
ground for each detector column. Wavelength calibration 
was done using the HeNeAr and the ThAr calibration spec- 
tra for the KPNO and NOT data, respectively, interpo- 
lating the wavelength solution between the nearest before 
and after calibration spectra. The average S /N of the spec- 
tra was 31 for KPNO and 25 for the NOT data, so while 
the KPNO telescope was much larger than the NOT, the 
longer NOT integrations, combined with a more efficient 
grating allowed similar S/N but higher spectral resolution 
in the NOT data. The spectra were normalized to the con- 
tinuum using the iraf task continuum which used a four- 
teenth order Legendre polynomial. The mean normalized 
KPNO and NOT spectra are shown in Fig. [2] on the same 
wavelength scale. The 10 and 25 s integration times lead to 
phase smearing that will reduce amplitudes of variation by 
0.85 and 5.2%, respectively, for a 140 s pulsation period. 
Additionally, our measurements are for disk-averaged ob- 
servations and intrinsic amplitudes may be higher locally 
within the unresolved surface. 

3. The frequency content of the pulsations 

Our selection of PG 1219 as a target is partially due to its 
year-to-year stability and this continued during our 2006 
observations. Table [3] provides the frequencies and ampli- 
tudes of pulsations detected during our four- year photomet- 
ric program and indicates that two additional frequencies 
were observed during 2004, but during single nights only. 
However, during our 2006 campaign, a new, low amplitude 
frequency persisted through both weeks of observations. 
These occasionally-observed frequencies are labeled as /5 
through /7 of Table [3] Pulsation spectra of photometry 
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Table 3. Frequencies and amplitudes detected in our annual photometric data sets. All frequencies are given in /iHz and 
amplitudes in mma with least-squares errors on the last digits in parentheses. 



ID 


2003 




2004 




2005 




2006 






Freq. 


Amp. 


Freq. 


Amp. 


Freq. 


Amp. 


Freq. 


Amp. 


n 


6721.484(10) 


2.8(1) 


6721.529(55) 


3.9(2) 


6721.481(158) 


2.2(4) 


6721.591(4) 


2.9(1) 


n 


6961.363(6) 


5.3(1) 


6961.355(33) 


6.6(2) 


6961.202(52) 


6.8(4) 


6961.489(2) 


6.4(1) 


/3 


7490.021(9) 


3.6(1) 


7489.727(38) 


5.7(2) 


7492.850(52) 


6.8(4) 


7490.833(2) 


4.5(1) 


/4 


7807.738(8) 


4.0(1) 


7807.771(29) 


7.4(2) 


7807.733(36) 


9.7(4) 


7807.857(2) 


6.1(1) 


/5 






7398.56(30) 


0.7(2) 










/6 














7744.6(45) 


2.1(4) 


/7 






8168.83(17) 


1.2(2) 












4200 4400 4600 

Wavelength (angstroms) 



Fig. 2. Average of all normalized KPNO (top) and NOT 
(bottom) spectra. 



which overlaps our spectroscopic data are shown in Fig. [3] 
Insets show the spectral window (FT of a single sine- wave 
sampled at the same times as the data) and smaller pan- 
els show the residuals after prewhitening by the four main 
frequencies. The solid horizontal (blue) line is the 4cr detec- 
tion calculated as four times the aver age value of the FT f or 
regions outside of areas of pulsation (jBreger et all 11993 ). 

Photometric amplitudes of fl through /4 show 20-40% 
variation over the course of four years, but during our spec- 
troscopic observations, the amplitudes remain essentially 
constant. The only significant change is /3, for which the 
amplitude increases by 14% in the week between the MDM 
and BO observations. From our photometric monitoring of 
PG 1219 we can conclude that there are four consistently 
detected frequencies and a fifth, low amplitude frequency. 

Since our cycle time and S/N were better for photom- 
etry, we could examine the photometry for nightly vari- 
ations in phase (defined as time of first maximum since 
JD=2453830.0 divided by the pulsation period). Such vari- 
ations can cause integrated data to indicate lower ampli- 
tudes or have multiple peaks which are not intrinsic to the 
star. The phases of the four main frequencies were consis- 
tent to within the errors of those provided in Table 4, and 
so we conclude that no phase variations occurred during 
the photometric or spectroscopic runs. However the phases 
do change between the MDM and BO observations. The 
phase of /3 changes by 13% while the others show marginal 
changes of 2, 4, and 3% for fl, /2, and /4, respectively. 




— w.^jiiit'ww 




1.5 

1 

0.5 



7000 7500 
Frequency (^Hz) 



Fig. 3. Pulsation spectra of photometric data simultaneous 
with the spectroscopic observations. Smaller panels have 
the four main frequencies removed and insets are the spec- 
tral windows plotted on the same frequency scale. The 
(blue) lines indicate the 4ct detection limit. 



3.1. Time-series spectroscopy 

The procedures we used a re similar to those described in 
iTelting fc Ostensenl (|2004l ) and summarized as follows: We 
calculated the radial velocity (RV) of the time-resolved 
spectra using the cross-correlation application FXCOR in 
I RAF. We produced nightly template spectra from the mean 
of individual spectra which FXCOR uses to do a Fourier 
cross-correlation between the template and individual spec- 
tra. We filtered the input to fit line features rather than 
large-scale trends in each spectrum, and a Gaussian was fit 
to the resulting cross-correlation function (CCF) of a best- 
fit size between 3 and 21 velocity bins around the maximum. 
We fit the CCF for the H/3 through H8 Balmer lines and 
the 4026 and 4471 A Hel lines. The radial velocity shifts 
are dominated by those of the strongest Balmer lines in 
the spectra. Figure Q] shows radial velocities obtained from 
one night each of KPNO and NOT data. The raw veloci- 
ties, shown in the top portion for each date, show jumps 
in velocity due to repositioning the target on the slit and 
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-50 
1! 



50 

50 
-100 - 



a Apr 



17 Apr. 



Time (Hours) 





7000 7500 
Frequency (^Hz) 



Fig. 4. Radial velocities obtained from cross-correlation of 
the time-resolved spectra. The top two panels show a run 
from KPNO both before and after fitting with polynomial 
segments and the bottom two panels show a run from the 
NOT before and after being fitted with polynomial seg- 
ments. All panels are plotted to the same relative scale. 



longer-term instrumental effects. These were corrected by 
fitting first or second-order polynomials to the segments 
between repositionings and are shown on the bottom por- 
tion for each date. All panels are plotted to the same scale 
but with varying velocity offsets. As these changes were of 
significantly longer time then the pulsations, there was no 
impact on pulsation velocities. 

Figure [5] shows the temporal spectra of the KPNO and 
NOT RV data separately. The solid (blue) line shows the 
4cr detection limit and the insets show the spectral window. 
Results from least-squares fitting for the amplitudes and 
phases are provided in Table|4]in §4. Apparent in the figure 
is the missing peak where /3 should be in the KPNO data. 
While the 4<r limit is nearly double that of the NOT data, 
/2, /3, and /4 have amplitudes in the NOT data which are 
sufficient to be detected in the KPNO data. /2 and /4 are 
both at least marginally detected, yet there is no trace of 
/3 at all. /I shows a noisy 'bump' in the appropriate place, 
but at just over 3cr, it is below the detection limit. 

We also compute the equivalent widths (EW) of the 
H/3 through H8 Balmer lines. For our EW measurements, 
we shifted the spectra using the cross-correlation veloc- 
ities to account for instrumental and pulsational veloc- 
ity shifts. The resulting standard deviation in velocity is 
0.07 km/s with an average CCF fit error of 5.34 km/s. To 
reduce the noise in the wings of the profiles, each point 
in the profile was weighted with its own depth with re- 
spect to the continuum, using the following non-standard 



Fig. 5. Pulsation spectra of radial velocity data. Insets are 
the spectral 

windows and (blue) lines indicate the Act detection limit. 



form, EW= \JY>{1- P{\)) x AA. No long-term trends 
were evident in the EW data, so no polynomial corrections 
were made. The results are shown in Fig. O fl is below 
the detection limit for both KPNO and NOT data, / 2 and 
/4 are easily recovered in both sets, and /3 is clearly de- 
tected in the NOT data, but just below 4cr in the KPNO 
data. As evident in Fig. [|5J NOT EW values were found to 
be systematically larger than for the KPNO data. Analysis 
using various resolutions and line-bin-widths has shown the 





KPNO I 




mind 








— 1 — ' — ' ' — 

NOT 


■\ ' ' 








Lyl 


Li 



Frequency (/xHz) 



Fig. 6. Temporal spectra of equivalent width data. 



amplitudes to be dependent on both the resolution and the 
number of pixels used to determine the line widths. Because 
of these dependencies it is not feasible to compare the am- 
plitudes between the two runs without some normalization. 
As the photometric and RV amplitudes of /4 are the clos- 
est between the runs, in Table [4] we have normalized the 
KPNO EW amplitudes to that of /4 from the NOT run, 
while we show their original amplitudes in Fig. [SJ 



4. Discussion 

4.1. EfFective temperature and gravity variations 

To measure changes in T r and log g for individual pulsation 
frequencies, we folded the spectra into 20 phase bins accord- 
ing to the pulsation periods. By folding over each pulsation 
period, we average out the variations caused by the other 
three periods, effectively isolating those of the folded pe- 
riod. The phase-binned spectra have S/N in the range 125 
to 160. We then fit these binn ed spectra with the LTE syn - 
thesis model atmospheres of iHeber. Reid. Wernerl ((2000). 
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The fi tting is a \ 2 process described by iBergeron et al.l 
(1 19921) and updat ed to include helium abundance fits by 
IS after et~a l. (1994). These model fits provide values of T e ff, 
logg, and helium abundance as a function of pulsation 
phase for each of the periods. As the helium abundance did 
not change during any of the pulsations, we fixed this to 
the mean value and refit the models, mildly decreasing the 
fit errors on T e g and log g. In Fig. [7] we illustrate the fitting 
procedure by showing the atmospheric fit to the template 
spectrum produced by combining all of the NOT spectra. 
The lines used for the fit are shown in the figure and include 
several Balmer lines and three He I lines. A single He II line 
is shown, but was not used in the fit. The fitted values for 
T c ff, log g, and N(He)/N(H) are also provided on the fig- 
ure and ar e nearly within the err o rs of those previously 
publi s hed (ICharpinet et ail. 120051: iHeber. Reid. Wernerl . 



120001 iKoen et all 119991 : lO'Toole fe HeberL 120061) . There 
are known systematic effects between LTE and NLTE 
model atmospheres w ith a nd without iron-group d iffusion 
(|Heber fe Edelmannl . 12004 lO'Toole fe HeberL I2006D and so 
we do not claim any increased accuracy in our values, but 
merely show the l ines t o illustrate the fitting procedure. 
lO'Toole fe Hebeil d2006h were able to overcome these in- 
consistencies using supersolar metal abundances, so we con- 
sider their measurements of T c r and log g as the most reli- 
able. In our case we are not really concerned with the actual 
values of T e g and log g but the changes that occur through 
a pulsation cycle. In this sense, minor differences from pre- 
viously reported values do not contribute to the fit errors 
in the variations. 



Hel 4472 
Hel 4922 
Hel 4026 

Hell 4686 (not fit) 
II 





T.„ = 33570 K (88 K) 
log g = 5.80 (0.013) 
N(He)/N(H) = -1.4673 (0.0323) 



-40 -20 20 

AA (Angstroms) 
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4.2. Fits of determined quantities 

To determine amplitudes and phases of the measurables (or 
model quantities), we used non- linear least-squares (NLLS) 
fits to measure the amplitudes and phases with the pulsa- 
tion frequencies fixed to the those derived from photometry 
(which match those in velocity). For photometry, RV, and 
EW, we fit the data themselves with timing corrected to 
the barycenter of the solar system. For T e g and log g, these 
are model fits to spectra folded over the frequency fixed 
to the photometric value. Then the variations in model T c g 
and log g were fit using our NLLS routine. The phase-folded 
data, along with the fits are shown in Figs. through [Tl] 
with the quantities provided in Table 2J Frequency-folded 
data were produced from the photometric and RV data 
for the plot, but these were not fitted for the quantities in 
Table H except as noted (/l and /3 RV data for KPNO 
and fl EW data for KPNO and NOT). Table H also pro- 
vides the 4(j detection limits from the data. No such values 
can be determined from the atmospheric model fits, though 
their formal errors have been folded into the NLLS formal 
errors shown in the figures and provided in the table. 

Since we know that the pulsations are phase-stable dur- 
ing our photometric observations, we can expect them to 
be phase-stable for the other measurables, and our results 
indicate it to be the case. Even the 13% change in phase for 
/3 between the MDM and BO runs is measured, to within 
the errors for RV, EW, T c $, and log g between the KPNO 
and NOT runs, even though some amplitudes show large 
variations. This indicates that our phases are measured re- 
liably. 



Fig. 7. Atmospheric models (dashed lines) fitted to the 
template spectrum of combined NOT data for several H 
and He I lines. The He II line is shown, but not used in the 
fit. Line identifiers and resultant atmospheric parameters 
are provided at the bottom. 



4.3. Amplitude ratios and phase differencing 

While we leave any attempts at mode identification to the 
second paper, which will properly simulate our results using 
pulsation models, it is still interesting to examine what we 
detected in such a context. It is expected that pulsations of 
various modes will behave differently between photometry, 
velocity, equivalent width, effective temperature, and grav- 
ity measurements, and that such differences (and similari- 
ties) will provide strong constraints on mode identifications. 
There has already been a model for this star matched to 
photometric data l |Charpinet et all 120051 ) which concluded 
that our four measured frequencies should be associated 
with three different pulsation degrees. As such, we should 
expect observable differences betw een our measured quan- 
tities. The mode identifications of ICharpinet et ail {2005) 
are provided in Fig. [T2l 

Figures [T2l through [T4l show amplitude ratios and phase 
differences of the observables using most of the available 
combinations. Figure [T^] shows those quantities directly ob- 
tained from the data, Fig. [T3] excludes photometry, but in- 
cludes spectroscopic quantities derived from models, and 
Fig U4l compares the photometry with the model quantities. 
Values for the NOT/BO data have solid errorbars while 
those for the KPNO/MDM data have dotted errorbars. We 
note that most of our quantities are within the la errors 



M. D. Reed et al.: Time-series Spectroscopy of PG 1219+534 



7 



Table 4. Results of photometric and spectroscopic variations in PG 1219. Notes: a indicates quantities were fit directly 
from the data, though they were below the 4ct limit. b indicates quantities that were fit from folded data (but not 
modeled). The KPNO EW amplitudes have been normalized to that of /4 of the NOT data. 



Measure 


Obs 


/I 


/2 


/3 


/4 


Limit 






Amplitudes 








Photometric 


MDM 


2.9(2) 


7.9(2) 


5.6(2) 


7.2(2) 


0.6 


(mma) 


BO 


3.3(2) 


7.8(2) 


8.0(2) 


7.1(2) 


1.1 


RV 


KPNO 


1.6 a (5) 


2.3(5) 


0.8° (0.5) 


2.2(5) 


2.2 


(km-s" 1 ) 


NOT 


1.4(3) 


3.1(3) 


3.0(3) 


2.4(3) 


1.2 


EW (1(T 3 ) 


KPNO 


13.3 a (5.4) 


40.0(5.4) 


26.2 a (5.4) 


39.4(5.4) 


28.2 




NOT 


16.2 a (3.8) 


34.7(3.8) 


34.8(3.8) 


39.4(3.8) 


18.2 


T c s 


KPNO 


46(26) 


175(44) 


58(30) 


159(44) 




(K) 


NOT 


59(26) 


122(34) 


130(34) 


111(35) 




logs 


KPNO 


12 (7) 


24(7) 


15(8) 


17(6) 




(l(r 3 cm-s- 2 ) 


NOT 


11(5) 


11(5) 
Phases 


11(5) 


16(5) 




Photometric 


MDM 


0.325(10) 


0.798(3) 


0.556(5) 


0.382(4) 






BO 


0.310(11) 


0.838(5) 


0.423(4) 


0.354(5) 




RV 


KPNO 


0.023 a (44) 


0.484(32) 


0.276 a (90) 


0.117(33) 






NOT 


0.060(33) 


0.545(15) 


0.142(15) 


0.066(19) 




EW 


KPNO 


0.837 a (65) 


0.310(22) 


0.074(33) 


0.860(22) 






NOT 


0.866 a (37) 


0.323(17) 


0.926(17) 


0.832(15) 






KPNO 


0.207(117) 


0.820(40) 


0.605(82) 


0.381(41) 






NOT 


0.316(58) 


0.850(44) 


0.441(41) 


0.349(52) 




logs 


KPNO 


0.392(113) 


0.794(51) 


0.444(72) 


0.384(69) 






NOT 


0.182(74) 


0.802(63) 


0.417(61) 


0.350(52) 





of all other quantities, when the largely instrumental EW 
offset between data sets is accounted for. While we do not 
attempt any mode identifications here, the similarities ex- 
hibited between the amplitude ratios suggest low degree 
modes with I < 2, which can behave similarly at most incli- 
nations. The I — 3 identification for /3 seems unlikely since 
this frequency behaves very much like the others. However, 
there are no published model amplitude ratios by which to 
judge our results. 

There are some expected relationships in phases that 
can be used to deduce the accuracy of the data as well. It 
is known that in low-degree modes, RV phases should be 
separated by one quarter (7r/2) from T c ff, logg, EW and 
photometric phases for adiabati c pulsations. This has been 
observed for other sdBV stars (jTelting fc Ostensenl [2004; 
iTillich et al.l.l2()()7i. and we see the same thing. It should 
also be expected that brightness, T e ff, and log g should all 
be in phase but EW should be in anti-phase to these ob- 
servables, and again, this is what we observe. Expected 
phase relations, deduced from atmospheric models are indi- 
cated as dashed lines in the figures. Our observations scat- 
ter tightly around these lines, indicating that our spectro- 
scopic reductions and model fitting are reliable. 



5. Conclusions 

During 2006 we collected simultaneous photometry and 
time-series spectroscopy from MDM and Kitt Peak 
Observatories and Baker Observatory and the NOT of the 
pulsating subdwarf B star PG 1219+534. In total, we ob- 
tained over 5000 spectra and nearly 15000 photometric 
measurements. We have recovered the four frequencies of 
the discovery paper (|Koen et al.l . ll999f ) in photometric, ve- 
locity, and equivalent width variations. 

From four years of photometry, we determine that 
PG 1219 has four consistent frequencies with occasional 




Fig. 8. Variations in T e g, log g, equivalent width, radial ve- 
locity, and photometry for fl — 6721 and /2 = 6961 fiHz 
of the NOT/BO data, along with non- linear least-squares 
fits (solid lines). Phases shown twice for clarity. 
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Fig. 9. Same as Fig. U for /3 
of the NOT/BO data. 



7489 and /4 = 7808 ^Hz Fig. 10. Same as Fig. [8] for KPNO/MDM data. 



low-amplitude transients. One such transient occurred dur- 
ing our spectroscopic runs, though we did not recover that 
frequency in any of our spectroscopically-determined quan- 
tities. The photometric amplitudes of /l, /2, and /4 remain 
stable over the course of our observations but the amplitude 
of /3 increased by 43% between our MDM and BO observa- 
tions. Similarly, the photometric phases were stable during 
our observations, but most varied between the runs. That 
of /3 showed the largest change in phase of 13% while /2 
and /4 had changes of 4 and 3%, respectively. Only /l had 
no phase changes, to within the errors. 

Our spectroscopic results were similar, in that only /3 
showed significant amplitude variations between the runs. 
It has a very low RV amplitude of 0.8 km/s in the KPNO 
data, which grows to an amplitude of 3.0 km/s in the NOT 
data. However, the differences in normalized EW ampli- 
tudes are within the errors for all frequencies. 

By folding the spectra over the pulsation period, we 
were able to fit atmospheric models to the higher S/N 
binned spectra to extract T e g and logg variations for all 
four frequencies. Again all pulsation amplitudes are con- 
sistent between the KPNO and NOT data, except for /3 
in T e ff which more than doubles. As such, /3 has signifi- 
cant amplitude increases of 43, 275, and 124% between the 
KPNO/MDM and NOT/BO observations for photometry, 
RV, and T c g , respectively. While there is a 33% increase in 
EW amplitude between the runs for /3, it is just within the 
la errors. 



We examined amplitude ratios and phase differences for 
various measurables for each frequency. We will interpret 
these more fully in a subsequent paper, but observation- 
ally, we can conclude that the ratios and differences do not 
differ greatly between frequencies. The simplest interpreta- 
tion of this result is that they are all low-degree (£ < 2) 
modes. We also compared phase differences to those ex- 
pected from adiabatic atmospheric models. Our measured 
differences match those expected, except for small depar- 
tures which are most likely caused by nonadiabatic effects. 

A large arsenal of quantities can now be applied to iden- 
tify the modes of pulsation. A future paper will match our 
measured quantities to those from perturbed synthetic pul- 
sation spectra to constrain the modes of each frequency. 
This work is also the first successful application of time- 
series spectroscopy to a "normal" amplitude sdBV star. 
All previous detections were of stars w i th unusually high 



photometric amplitudes (10 'Toole et all. l2000l [200 



2005: IWoolf et all l2002t 



Tclting & Ostcnscn, 200 



2003 



2006 



iTillich et all l2007f ) and so our work indicates that such 
studies can provide useful measurements for the majority 
of sdBV stars. 
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Fig. 11. Same as Fig. U for KPNO/MDM data. 
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Fig. 12. Amplitude ratios versus phase differences for fitted 
observables. Solid errorbars indicate NOT/BO data while 
dotted errorbars indicate KPNO/MDM data. Dashed lines 
indicate expected phases fr om model simula t ions a nd mode 
identifications are those of ICharpinet et all ({2005) . 
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Fig. 13. Same as Fig. Q2] for quantities determined solely 
from spectroscopy. 
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Fig. 14. Same as Fig. [T2lcomparing photometry and model- 
fitted quantities. 
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